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Abstract Recent experiments have demonstrated a remarkable progress in im- 
plementing and use of the Particle Image Velocimetry (PIV) and particle tracking 
techniques for the study of turbulence in 4 He. However, an interpretation of the 
experimental data in the superfluid phase requires understanding how the motion 
of tracer particles is affected by the two components, the viscous normal fluid and 
the inviscid superfluid. Of a particular importance is the problem of particle inter- 
actions with quantized vortex lines which may not only strongly affect the particle 
motion, but, under certain conditions, may even trap particles on quantized vortex 
cores. The article reviews recent theoretical, numerical, and experimental results 
in this rapidly developing area of research, putting critically together recent re- 
sults, and solving apparent inconsistencies. Also discussed is a closely related 
technique of detection of quantized vortices negative ion bubbles in 4 He. 
PACS numbers: 67.40.Vs,47.37.+q,47.27.-i 



1 Introduction and plan of the review 

In this review we will be mostly concerned with the motion of small solid particles 
in turbulent 4 He. This new and rapidly developing area of research has been ini- 
tiated by the recent success of implementation of the Particle Image Velocimetry 
(PIV) and the particle tracking techniques in superfluid helium, see Donnelly et 
al.—, Van Sciver and co-workers 2 - 3 A 5 i 6 ; and publications of the group involving 
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Paoletti, Bewley, Lathrop, Sreenivasan and their co-workers 7 i 8 i 9 i 10 i 1 1 . (The only 
difference between these two, otherwise identical, techniques is that the results 
of PIV yield the local average particle velocities obtained by calculating cross- 
correlations of particle ensembles and, therefore, result in the fluid-like, smooth 
velocity field, while the particle tracking aims at investigation of individual parti- 
cle trajectories.) These are perhaps the only two well developed techniques which 
can identify the flow patterns in turbulent superfluid helium. The PIV has been a 
standard technique in classical fluid dynamics for several decades (see, for exam- 
ple, the book by Raffel et al. 12 ). In experimental studies of classical turbulence, 
one can be confident that the motion of sufficiently small particles will reveal the 
details of turbulent motion of the viscous fluid. In superfluids, the two-fluid nature 
of 4 He makes an interpretation of PIV measurements much more difficult: A solid 
particle interacts with both the normal fluid and the superfluid; moreover, the parti- 
cle may interact strongly with quantized vortex lines in the superfluid component, 
and even become trapped on them. 

The mechanism of trapping of solid particles is essentially the same as that of 
trapping the charge carriers on quantized vortices in 4 He, primarily because e.g. 
the negative ion (electron) forms around itself an almost macroscopic bubble of 
radius 12 — 20 A from which helium atoms are excluded. The ion trapping tech- 
nique was used for detection of quantized vortices since late 50s, and, although 
not being suitable for studying the normal fluid patterns in 4 He, can be consid- 
ered as the technique closely related to the PIV and the particle tracking methods. 
We will review the development of the ion trapping technique in Section 7 of this 
article. 

The aims of experimental physicists working in the area of superfluid turbu- 
lence are to understand the vortex tangle (quantum turbulence), to measure veloc- 
ity fluctuations in both the normal fluid and the superfluid, and possibly to make 
comparisons between quantum turbulence and ordinary turbulence. Until recently 
the major difficulty was lack of direct flow visualization technique near absolute 
zero. Existing methods, such as measurements of temperature differences to de- 
tect extra dissipation, ion trapping, measurements of second sound attenuation by 
quantized vortices, etc. probed only the vortex line density L (total vortex length 
per unit volume) averaged over an experimental cell and did not reveal turbulent 
velocity fluctuations (although a remarkable resolution has been recently achieved 
in second sound measurements by Roche et al— of the local tangle density). 

Experimental breakthrough was made in 2002 when the PIV technique was 
successfully introduced in 4 He (see references in the first paragraph of this Sec- 
tion). This technique, which has been standard in classical fluid dynamics for 
many years, is based on injecting many small particles into the liquid. Two im- 
ages are produced using short laser pulses of different frequencies (corresponding 
e.g. to the green and red colors) focused into a narrow sheet and separated in time 
by a few milliseconds. Software then analyzes the images and identifies green and 
red dots corresponding to the same particle at the two different times. In this way 
the observed distance between the corresponding dots yields the component of the 
local velocity in the plane of the light sheet. 

In classical fluids, provided the particle size is sufficiently small (in turbulence 
studies, "small" means smaller than the Kolmogorov length), the dominating force 
acting on the particle is the viscous drag force, so that a researcher can be con- 
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fident that small particles trace the fluid motion (in particular, turbulent velocity 
fluctuations) rather well. 

Because the viscosity of liquid helium is very low the Kolmogorov scale in 
the turbulent normal fluid can be very small, so that it is essential that the particles 
used in visualization experiments are as small as possible. In the cited experimen- 
tal works the typical particle diameter was of the order of 1 fim. 

What do tracer particles trace in superfluid helium? One may expect that, al- 
though the viscosity of 4 He is low, the dominating force exerted on the particle 
by the fluid will still be the viscous drag, so that solid particles should trace the 
normal fluid. However, this is not always true: due to the two-fluid nature of su- 
perfluid helium, the particles interact not only with the normal fluid, but also with 
the superfluid component through inertial and added mass forces; moreover, the 
particles interact strongly with quantized vortices in superfluid and may even be 
trapped on superfluid vortex lines. Therefore, if we want to interpret results of 
PIV and particle tracking measurements correctly, we must answer first the ques- 
tion asked in the beginning of this paragraph. 

This article is divided in sections where theoretical and numerical models of 
increasing complexity are compared to each other and to experimental results. The 
first model, described in Section 2, is the one-way coupling model. In Sections 2. 1 
and 2.2 we derive the governing equations of motion of particle tracers in the pres- 
ence of two imposed fluids, the viscous normal fluid and the inviscid superfluid, 
under the assumptions that the particles do not disturb the flow, are smaller than 
any flow scale of interest, and do not become trapped in vortices. The one-way 
coupling model allows us to discuss the problem of the stability of particles' tra- 
jectories, which is relevant to the visualization a pure superflow. The one-way 
coupling model is powerful enough to derive the general principles which lead 
to particles being trapped on vortices. In Section 2 we present the experimental 
evidence for this trapping (Section 2.3), and show how the mutual friction affects 
the motion of particles near vortex cores (Section 2.4). Three case studies are 
discussed of particle trajectories near vortices: vortex ring, thermal counterflow 
tangles, and vortex tangles at low temperatures (Section 2.5). 

Section 3 introduces the more sophisticated (and computationally more expen- 
sive) two-way coupling model; in this model the back reaction of the flow on the 
particle is taken into account, and the dynamics of the particle-vortex interaction 
and the trapping can be studied in great detail, including what happens at the sur- 
face of the particle (which we assume to be spherical for simplicity). Section 3.1 
contains the mathematical formulation of the two-way coupling model; Sections 
3.2 and 3.3 are devoted to the numerical calculations of typical vortex-particle 
interactions. 

Section 4 makes use of the results of the numerical simulations described in 
Section 3 to derive a simpler analytical model which explains in a quantitive way 
experiments performed in Florid a 2 i 3 i 4 i 5 i in which heavy particles fell through a 
tangle of vortices. The two-way coupling model is also applied to particles moving 
in tangles generated by a thermal counterflow: these numerical results are applied 
to particle-tracking experiments performed in Maryland 7 ' 8 ^ 10 ' 11 . By considering 
numerical calculations at small and high values of the vortex line density, we solve 
the apparent disagreement between the Florida and Maryland experiments: the 
discussion will reveal that the two experiments refer to two different regimes, 
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which are both explained by the two-way coupling model. In Section 4 we also 
discuss the experimental observation of flow structures observed behind and in 
front a cylindrical obstacle set in the middle of a counterflow channel, and propose 
a simple analytical two-dimensional model which accounts, at least in principle, 
for the qualitative existence of these flow structures. 

Section 5 describes the most recent experimental discoveries obtained using 
tracer particles: the observation of turbulent boundary layer flows, the visualiza- 
tion of individual vortex reconnections, and the measurement of velocity statis- 
tics. Section 7 reviews other techniques based on trapping ions and imaging He2 
molecules; these techniques share important principles (but not the size of the 
trapped object) with PIV and particle tracking techniques. Section 8 contains the 
final discussion. 



2 One-way coupling model of particle motion in turbulent 4 He 

We begin answering this question with a relatively simple, "one-way coupling" 
model which follows the approach standard in classical two-phase turbulence 
studies. We will generalize to the two-fluid model of superfluid helium the equa- 
tions of motion of a solid spherical particle of radius a p in a nonuniform flow of 
classical fluid. We assume that 1) the presence of particles does not modify the 
turbulence, 2) the flow velocities vary little over a distance of the order of particle 
size, and 3) particles do not interact strongly with quantized vortex lines and are 
certainly not trapped on these lines. These assumptions require that the particle 
size be much smaller than both the Kolmogorov lengthscale, b-q in the normal 
fluid, and the mean intervortex distance, £ = L -1 / 2 in the superfluid. Below we 
formulate, under these assumptions, the equations of particle motion. 



2. 1 Fluid-particle interaction 

We start with the fluid-particle interaction in classical liquids. We consider a spher- 
ical solid particle of radius a p in the nonuniform flow. Let the ambient (that is, in 
the absence of the particle) fluid velocity field be v(r, t). According to the as- 
sumptions formulated above, the size of the particle is much smaller than the flow 
lengthscale, A, i.e. a p <C A, so that we can introduce a small parameter 

£ = a / ,||Vv||/|v-u p | < 1, (1) 

where u p is the particle velocity. 

2.1.1 Fluid-particle interaction in the inviscid nonuniform flow 

As shown by Auton et al. 14 , the total force acting on the particle in the nonuni- 
form, inviscid flow can, under assumption (flj, be represented in the form 

F = + + F^ , (2) 
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and C is the added mass coefficient (for spherical particle C 




(5) 



where ft) = V x v is the vorticity, represents the lift force arising due to stretching 
of vortex lines in the vicinity of the sphere's surface. 



2.7.2 Fluid-particle interaction in the nonuniform viscous flow 

We will consider the motion of a solid particle in the viscous fluid assuming that 
the particle Reynolds number is small: 



where V is the kinematic viscosity. Note that small particle Reynolds numbers 
are typical of PIV and particle tracking experiments, both in classical fluids and 
in superfluid helium. We will consider the particle motion under the assumptions 
formulated above (in particular that the parameter £ introduced by formula ([]]) is 
small). 

Detailed analysis of the forces acting on the particle in the nonuniform viscous 
flow can be found in works by Maxey and Riley—, Mei— , and Kim, Elghobashi, 
and Sirignano 17 . The total force acting on the particle can be represented as a sum 
of several contributions, i.e. the gravity, viscous drag, the inertial and added mass 
force, Faxen correction arising due to the local non-uniformity of the ambient 
flow, the Saffman lift force arising due to the local shear, and the Magnus lift 
force arising due to rotation of the particle. 

The main contribution, dominating in most particulate flows, is the viscous 
Stokes drag force: 



Surprisingly, in the viscous flow the inertial and added mass forces, and 
F( a ) are determined by the same formulae ©-(O as for the inviscid flow, with the 
same added mass coefficient, C = \ for the spherical particle. 



Re p = a p \\-u p \/v < 1 , 



(6) 



F (a) =6na p pv{\-Up). 



(7) 
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2.1.3 Fluid-particle interaction in 4 He 

To determine the force acting on the particle in 4 He, we simply add together all 
the forces exerted by the normal fluid and the superfluid. Since the superflow is 

potential, F^ m ' is identically zero, provided the particle does not become trapped 
on quantized vortex lines. As shown by Poole et al.—, for the flow properties and 
physical parameters typical of the normal component of 4 He, the Faxen correction, 
the history and lift forces can be neglected provided the particle Reynolds number 
and the parameter £ are small and the particle size is significantly smaller than 
the Kolmogorov length. Therefore, the total force acting on the particle can be 
approximated as 

F = F® +Fi d) +F« +Fi a) +F? +F< a) , (8) 

where the subscripts n and s refer to the normal fluid and the superfluid, respec- 
tively. Here F® is the combination of gravity and buoyancy, and F^ is the vis- 
cous drag force exerted by the normal fluid: 

V®=#(p p -p)g, Fi d) =6^M«(v„-u p ), (9) 

where ji„ is the viscosity of 4 He. Substantial derivatives required for determining 
the inertial and added mass forces in the normal and superfluid components are 
now defined by formula © with v replaced by v„ and \ s , respectively. 



2.2 Lagrangian equations of particle motion 

We arrive at the following equation of particle motion—: 

Pp$ ^ = 67ta p n n (v n - Up) + $(p p - p)g 

This equation must be considered together with the kinematic equation 

dr/dt=u p , (11) 

where r = r(f) = (x(t),y(t), z(t)) (= r p (t)) should be regarded as a Lagrangian 
trajectory of the solid particle. Eqs. d 1 Ob and dTTT) constitute the closed system for 
the unknown particle position and velocity, r(?) and u p (t), respectively. 

Often (although not always) in the PIV and particle tracking experiments neu- 
trally buoyant particles (with p p = p) are used in order to eliminate unwanted 
effects of gravity on the particle motion. For neutrally buoyant particles, Eq. dlOl) 
can be written in a more concise form 

dw.nl, , 3 / D\ n D\ s \ 

-dr = r^-^ + 2 P - { p "^r +ps ^)^ (12) 
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where 

po=P P + p/2 = 3p/2 and t = 2a 2 p p / (9pL n ) . (13) 

The parameter T, which shows how quickly the particle adjusts its motion to the 
viscous flow, plays an important role in the study of particle motion in fluids, and 
is commonly known as either the particle response time, or viscous relaxation 
time. 

Analyzing Eqs. dTTb and dl2l l, Poole et al. 18 showed that, provided t/t/ <C 1, 
where ty is the timescale of the fluid motion (e.g. the Kolmogorov time), the 
neutrally buoyant particle tracks the motion of the normal fluid. If t/tf S> 1, 
the particle moves with a velocity corresponding to the total current density, j = 
P«v„ + p s v s . 

However, two important issues may invalidate these conclusions: 1) instability 
of particle trajectories, and 2) trapping of particles on superfluid vortex lines. The 
analysis of particle trapping will require more elaborate, self-consistent, two-way 
coupling model which would account for deformation of the vortex filament by 
the approaching particle, including possible reconnection of the vortex with the 
particle surface (see below Section [3 . 



2.2.7 Instability of particle trajectories 

Instability of Lagrangian trajectories of the neutrally buoyant particle in the clas- 
sical viscous fluid was discovered and studied relatively recently by Babiano et 
al.—. To illustrate such an instability in turbulent superfluid at finite temperature 
such that both the normal fluid and the superfluid component are present, we con- 
sider the simplest case assuming that v„ « v, =: v down to the length scales com- 
parable to the vortex line spacing. The Lagrangian equations of particle motion, 
which become 

+ — ^-u (14) 
dt ~ % Dt 1 dt Up > (14) 

have, provided u p (0) = v(0) and r p (0) = iy(0), a formal solution u p (t) = v(ry, f), 
r p (t) = ry(f), where iy(f) is a trajectory of the fluid point. Therefore, it would 
seem natural to conclude that neutrally buoyant particles follow the fluid exactly. 

However, let us consider now the particle motion in the so-called ABC (Arnold- 
Beltrami-Childress) flow, frequently used as the simplest model of turbulence, 
whose velocity field is 

v x = Asin(27rz) + Ccos(2?ry) , v y =Bsm(2nx) +Acos(27Tz) , 

v z = Csin(27Ty) +Bcos(27D:) . (15) 

The time sequence 18 illustrating the position of large number of neutrally buoy- 
ant tracer particles in the ABC flow, starting from their uniform distribution, is 
shown in Fig. [TJ It can be seen that particles do not follow the fluid, but instead, 
due to instabilities of their trajectories, travel from the regions of high vorticity to 
the regions of high rate of strain. (In classical multiphase fluid dynamics similar 
phenomenon for particles heavier than the fluid has been studied in detail. How- 
ever, the mechanism of segregation of heavy particles is quite different from that 
of neutrally buoyant particles.) Note, though, that this instability develops rather 
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Fig. 1 Positions of tracer particles at different timesia. From Poole, Barenghi, Sergeev, and 
Vinen, Phys. Rev. B, 71, 064514, (2005). Reprinted by permission, ©2005 American Physical 
Society. 



slowly: the last frame of Fig. Q] recovering the lines of minimum vorticity and 
maximum rate of strain, corresponds to 3 times of turnover of the ABC flow. 

Instability of particle trajectories is more pronounced in the case T — > when 
the normal fluid is absent. The equations of motion of neutrally buoyant particle 
become 

-^ = ^ + (v i .V)v s = -V^ ) -±=u p> (16) 

and have a formal solution U„(t) = \ s (r s (t), t), r p (t) = r,(f), where dr p /dt = 
\ s (r s ,t), so that, in the case of very low temperature when the normal fluid is 
absent, it can be expected that neutrally buoyant particles trace the superfluid. 
However, this is not the case either. Consider the simplest case of the neutrally 
buoyant particle moving around a single, stationary, straight vortex line. In cylin- 
drical polar coordinates (r p , 9 p ) the equations of particle motion are: 

f p - r p (0 2 p = - k 2 /{In 2 ?!,) , 2(t)r p + r p d> p = , (17) 

where (O p = 6 p . If at the initial moment the particle velocity does not coincide 
exactly with the fluid velocity, the particle will spiral either outwards, or inwards. 



9 



30 




-3C 1 ' ! ' 1 ' 1 

-30 -20 -10 10 20 30 



Fig. 2 Trajectories- of neutrally buoyant solid particle (dashed line A to C) and superfluid 
particle (solid line A to B) around three vortices moving along the closed orbit shown. From 
Sergeev, Barenghi, Kivotides, and Vinen, Phys. Rev. B, 73, 052502, (2006). Reprinted by per- 
mission, ©2006 American Physical Society. 



We arrive at the conclusion which remains valid in the general case of particle mo- 
tion at temperature T —* (Sergeev et al.—): unless the initial velocity of neutrally 
buoyant particle matches exactly the velocity of the fluid point, the trajectory of 
the solid particle deviates significantly from the trajectory of the fluid point. Such 
an instability is amplified by any macroscopic mismatch between the velocity of 
the superfluid and the velocity of the particle at the beginning of the experiment. 
This is further illustrated by Fig. |2] showing the trajectories of solid and fluid par- 
ticles around three vortices. 



2.3 Trapping of particles on quantized vortices: mechanism and experimental 
evidence 

2.3.1 Mechanism 

Why can the superfluid vortex trap the particle? To answer this question, we have 
to take into account a possibility of reconnection of the quantized vortex to the 
surface of moving particle (a more detailed analysis of the mechanism of this 
process will be given below in Section [3]). 

Below three different versions, or, rather, three different ways of explaining 
the reason for particle trapping are suggested. 

1°. Imagine that the quantized vortex reconnects symmetrically to the surface of 
spherical particle, as shown in the right part of Fig. [3] The kinetic energy of the 
flow field created by the straight quantized vortex can be easily calculated, and in 
the symmetric configuration is reduced by 

AE= P -4^Xn^, (18) 

where t; w 10~ 8 cm is the vortex core radius. Note that assuming a p » ^ this re- 
sult follows from the substitution energy calculated by Parks and Donnelly — for 
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Fig. 3 Asymmetric reconnection of the vortex to the particle surface creates a force restoring the 
symmetric particle-vortex configuration--. From Sergeev, Barenghi, and Kivotides, Phys. Rev. 
B, 74, 184506, (2006). Reprinted by permission, ©2006 American Physical Society. 

the ion bubble trapped by the quantized vortex line, see below Sec. Q and for- 
mula (138b therein. Formula (118b determines the kinetic energy which the particle 
would require to break free from the symmetric vortex configuration shown in 
Fig.EKright). 

2°. The flow field of the vortex creates a pressure gradient 

Vp = -p,(v,.V)v,= ^v(i) (19) 
attracting the particle to the vortex. 

3°. If the particle-vortex configuration is symmetric, as in the right part of Fig. [3] 
then, obviously, the force acting on the particle is zero. Now imagine that this 
symmetry is perturbed as shown in the left part of this figure. The superfiuid com- 
ponent is inviscid, so that the following arguments apply based on Bernoulli's 
integral: the fluid velocity on the right side of the particle surface, where the two 
vortex strands are closer, is greater than that on the left side, and, therefore, the 
pressure is bigger on the left side of the sphere. This provides a net force restoring 
the symmetric particle-vortex configuration. 

2.3.2 Evidence of particle trapping 

Experimental evidence comes from the publication of Bewley, Lathrop and Sreeni- 
vasan 7 (see also more recent paper s W 1 ). Fig. [4] published in the cited paper 7 , 
shows that the researchers, using the PIV technique, actually "painted" quantized 
vortices by tracer particles. 

Another evidence comes not from experiments, but from the quantum calcu- 
lation based on the Gross-Pitaevskii equation. Berloff and Roberts 23 studied an 
interaction between the negative ion and the quantized vortex. In superfiuid he- 
lium, the negative ion forms around itself a bubble of diameter approximately 
16 x 10~ 8 cm (an order of magnitude higher than the size of the vortex core, 
t, « 10~ 8 cm) and, therefore, can be treated, for our purpose, as a particle (albeit 
very small). Results of Berloff and Roberts' calculation are illustrated in Fig. [5] 
the ion bubble approaches the vortex, which deforms ((b) and (c)) and then recon- 
nects to the particle surface (d); the reconnection excites Kelvin waves ((f), (g), 
and (h)) which carry away the energy, and, eventually, the bubble-vortex configu- 
ration relaxes and the ion becomes trapped on the vortex core (i). 
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Fig. 4 PIV visualization 7 by Bewley, Lathrop, and Sreenivasan (Nature, 441, 588, (2006)) of 
quantized vortex cores: (a) - above A transition; (b) and (c) - branching filaments tenth of mK 
below X transition; (d) - regrouping along vertical lines in steady rotating apparatus. Reprinted 
by permission, ©2006 Macmillan Publishers Ltd. 

2.4 Mutual friction and trapping 

At this point, the question can be asked whether it is worth or not further exploit- 
ing the one-way coupling model which neglects any influence of the particle on 
the vortex evolution. The answer is "yes": there are cases where trapping events 
are not very frequent, so that a useful information about the behaviour of tracer 
particles can be obtained by ignoring their trapping on vortex cores. We will also 
show that the mutual friction between quantized vortices and the normal fluid can 
prevent trapping. 

In the vicinity of the vortex core, the mutual friction induces, in the normal 
fluid, the vortex dipole whose typical lengthscale is expected to be about 100 jim. 
The results of numerical calculation by Idowu et al. 24 of the dipole-like normal 
fluid disturbance are shown in Fig. [6] This normal flow disturbance can deflect the 
tracer particle which otherwise would have collided with, and possibly trapped 
by the vortex. Typical trajectory, calculated by Sergeev et al. 25 , of the particle 
moving from right to left and interacting with the superfluid vortex and normal 
fluid disturbance is shown in Fig.|7]by the solid line. For comparison, the trajectory 
calculated without taking into account the normal fluid disturbance is shown by 
the dashed line. This trajectory leads to the collision with the vortex core located 
at the origin. 

Normal fluid vortical structures induced by the mutual friction were predicted 
by Hall and Vinen 26 and Kivotides, Barenghi, and Samuels 27 , but so far, because 
of low resolution of experimental techniques, there was no direct experimental 
proof of existence of the normal fluid disturbances. Owing to the much higher 
resolution, the PIV and particle tracking techniques can provide such an evidence. 
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Fig. 6 Velocity field of the normal fluid due to the mutual friction forcing of a single superfluid 
vortex—. (A): velocity arrow plot; (B): streamlines. From Idowu, Willis, Barenghi, and Samuels, 
Phys. Rev. B, 62, 3409, (2000). By permission, ©2000 American Physical Society. 
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Fig. 7 Particle trajectories in the presence (solid lines) and absence (dashed line) of the normal 
fluid disturbances. From Sergeev, Wang, Meneguz, and Barenghi, J. Low Temp. Phys., 146, 417, 
(2007). Reprinted by permission, ©2007 Springer. 



Fig. 8 Observation by Bewley, Lathrop, Sreenivasan, and Paoletti of particle trajectories pos- 
sibly indicating an influence of normal fluid disturbances induced by the mutual friction. 



Perhaps the first experimental confirmation of existence of normal flow structures 
induced by the mutual friction has already been found. Fig. [8] shows a typical 
particle trajectory observed by Bewley, Lathrop, Sreenivasan, and Paoletti — , by 
means of the particle tracking technique, in the thermal counterflow. The particle 
trajectory has characteristic 'kinks' resembling the 'deflected' trajectory shown 
by the solid line in Fig. [7] 

In the remaining part of this Section we will analyze, neglecting trapping, three 
examples of particle motion in 4 He. 
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Fig. 9 Vortex ring and solid particles' configurations before and after the ring has passed through 
the particulate sheet 29 . Adapted from Kivotides, Barenghi, and Sergeev, Phys. Rev. Lett., 95, 
215302, (2005). By permission, ©2005 American Physical Society. 

2.5 Case studies 

2.5.1 Vortex ring propagating against a particulate sheet 

This study 29 led to a proposal of a simple experiment (not performed yet) which 
could, in principle, justify the use of PIV technique for measuring instantaneous 
normal fluid velocity patterns. We consider a single vortex ring propagating nor- 
mally to a plane sheet of neutrally buoyant particles, see Fig.[9](left). 

In the considered model, the mutual friction between the normal fluid and the 
superfluid vortex is taken into account leading to emergence of the normal fluid 
disturbances described earlier. The motion of the superfluid vortex whose core is 
defined parametrically by X(s, t) is governed by the following equation derived 
by Idowu et al.—. 

dX/dt = \ e = h\ s +h,X' x (v„ - V s ) - /z**X' x(X'xv„), (20) 

where X' = dX/ds, h(T), h*{T) and h t>t (T) are the known mutual friction coef- 
ficients, and the vortex-induced superfluid velocity V s is given by the Biot-Savart 
integral 



The motion of incompressible (V • v„ = 0) normal fluid is governed by the equation 



where F is the mutual friction force per unit volume. This force is determined as 
the sum of the drag force and the Iordanskii force, f, per unit length: 




(21) 



^ + (v„ • V)v n = -- Vp + vV 2 v„ + -F , 



P P 



(22) 



i = p s K[d^X'x (X' x (v„-V s ))-^X'x (v„-V s )]. (23) 
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Fig. 10 Histograms of the cosine of the angle between v„ and u p (left) and of the magnitude 
of relative velocity, [v„ — u p | (right). The vertical axes are divided by 1000. From Kivotides, 
Barenghi, and Sergeev, Phys. Rev. Lett., 95, 215302, (2005). Reprinted by permission, ©2005 
American Physical Society. 



Here the new mutual friction coefficients, d*(T) and (T) are known and can be 
expressed explicitly through h(T), h*(T) and h„(T). 

Trajectories of neutrally buoyant particles were found by numerical integra- 
tion of Eqs. Jill) and (112b . Configurations of solid particles before and after the 
ring has passed the particulate sheet are shown in Fig.|9](right). Fig.[l0]shows the 
histograms of the angle and of the relative difference between v„ and u p ; it can 
be seen that v„ and u p , to a very good degree of accuracy, are identical both in 
magnitude and direction. Moreover, the calculation showed that trapping events 
are relatively rare - only 42 out of 900 particles approached the ring to a dis- 
tance smaller than three particle diameters. This enables us to conclude that, in 
the proposed experiment, the measurement of particle velocities can provide di- 
rect information about instantaneous normal flow patterns. 

Later, based on the self-consistent, two-way coupling model, Kivotides and 
Wilkin 31 performed more elaborate study of interactions between solid particles 
and vortex rings, see below Sec. 13.41 

2.5.2 Particle motion in thermal counter/low 

In this example 3 — we will be concerned with the T-I state of 4 He turbulence such 
that the vortex tangle in the superfluid component is present but the normal flow 
is laminar. For simplicity, the normal flow is assumed uniform, v„ = const. The 
superfluid velocity can be represented as \ s = \f + V s , where the mean (coun- 
terflow) superfluid velocity, is linked with the normal fluid velocity by the 
relation pvf + p„v„ = 0, and V s is the fluctuating superfluid velocity induced 
by the vortex tangle. The dynamic vortex tangle is modeled, in the periodic box, 
taking into account the mutual friction between the normal fluid and quantized 
vortices, as well as the Biot-Savart interaction between vortex filaments. An in- 
fluence of superfluid vortices on the motion of normal fluid is neglected. Based 
on Eqs. (1111) and (112b . the motion of neutrally buoyant tracer particles of diam- 
eter dp = 6.25 x 10 _3 cm was calculated after the vortex tangle has reached the 
statistically steady state. Calculations, performed at temperature T = 1.3K for 
|v„| = 1.1417 and 0.6058 cm/s (corresponding values of the counterflow heat flux 
are q= 1.07 x 10~ 3 and 4.57 x 10~ 4 J/(cm 2 • s)), and at temperature T = 2.171 K 
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Fig. 11 Histograms 32 of cos(V s , u p )\. From Kivotides, Barenghi, and Sergeev, Europhys. 
Lett, 73, 733, (2006). Reprinted by permission, ©2006 EDP Sciences. 



for |v„| = 0.011 83 cm/s (q = 0.125 J/(cm 2 • s)), showed that the particle velocity 
is very narrowly peaked around the constant normal velocity. 

Three histograms of Fig.[TT] corresponding to the three cases considered above, 
show the absence of alignment (left) between particle velocities and the superfluid 
velocity, V s induced by the vortex tangle, and illustrate anisotropy in V s (center 
and right). 

The above two examples enable us to expect that, in the range of parameters 
typical of PIV measurements, neutrally buoyant particles should trace the normal 
fluid well, and, provided strong interactions of particles with quantized vortices 
can be neglected, particle velocity fluctuations induced by interactions between 
particles and quantized vortices should be relatively small. However, these con- 
clusion can be invalidated by trapping of solid particles on quantized vortex cores. 
This phenomenon will be addressed below in the Sec. [3] 



2.5.3 Particle motion in a vortex tangle at very low temperature 

Below the results discussed in Section [331 will suggest that at temperature T < 
0.5 K, when the normal fluid is practically absent and the damping force on the 
particle can be neglected, trapping of neutrally buoyant solid particles on quan- 
tized vortices can, most likely, be ignored and the motion of solid particle can be 
modeled by a simpler, one-way coupling model. However, at these temperatures 
the presence and motion of the particles still affects the vortex filaments, so that a 
certain modification should be necessary of the one-way coupling model. In order 
to understand some features of the particle motion in the vortex tangle at such a 
low temperature, we start with a simple, two-dimensional model of the tangle^. 
In such a model the vortex lines become vortex points, and the Biot-Savart law 
reduces to a simple statement that each vortex point moves as a fluid point in a 
flow field created by all other vortices. (Such a two-dimensional system of vortex 
points is known as the Onsager's point vortex gas^i.) 

We consider, in the periodic box, the motion of neutrally buoyant particles in 
the system of vortex points of random polarity. We neglect trapping of particles 
on quantized vortices as well as any influence of particles on the motion of vortex 
points. 

In this approximation the particle motion is governed by Eqs. ( 116b . and the 
following problem can be immediately identified: in the case where the posi- 
tion of the particle coincides with that of the vortex point, the pressure gradient 
force, —Vp becomes unphysically singular. Based on the mechanism, described 
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Fig. 12 Left: trajectories of a solid particle (dashed line) and of a fluid point (solid line) in a 
system of 20 vortex points. Right: projection of the particle trajectory in the three-dimensional 
vortex tangle 33 . From Kivotides, Sergeev, and Barenghi, Phys. Fluids, 20, 055105, (2008). 
Reprinted by permission, ©2008 American Institute of Physics. 

below in Sec.[3]of the particle-vortex collision at very low temperature (in the ab- 
sence of damping force), we will resolve this manifest difficulty by modifying the 
model d 1 6b as follows. 

The cause of the problem is that in the real, three-dimensional tangle at temper- 
ature T < 0.5 K, even when the particle breaks through the vortex, it nevertheless 
reconnects with the vortex filament when the distance between the particle and the 
vortex core becomes of the order of particle radius. Since the vortex line attached 
to the particle is necessarily orthogonal to the particle surface, the reconnection 
results in a dramatic decrease of the force exerted on the particle; this force is 
zero when the particle- vortex configuration is symmetric. Then, the one-way cou- 
pling model can be modified by assuming that there exists a force-free region for 
r p < a c , where a c ~ 0(a p ) is a cut-off distance. 

In the two-dimensional calculation we set ac = a p . Trajectories of an inertial 
particle and a fluid point are illustrated in Fig. Q~2] (left) for a system of 20 vortex 
points set, initially, at random locations. The solid particle starts its motion at the 
point A where it has the velocity equal that of the fluid particle, but very quickly 
the trajectory of the solid particle looses any resemblance to the trajectory of the 
fluid point; moreover, the trajectory of the solid particle soon acquires the ballistic 
character. The reason for such a behaviour of the solid particle is the instability of 
its trajectory as discussed in Sec. 12.2.11 

The evolution of particle velocity with time reveals another, rather unexpected 
feature: although |u p (0)| = |v s (0)|, the magnitude of particle velocity quickly in- 
creases above that of the fluid point, |v s (f)| and eventually saturates remaining 
larger than the average value of |v s | at all times. The average saturated particle 
velocity satisfies the scaling 



(u p )~(v s }^/I/a~ c , (24) 

where I is the intervortex spacing, and (v s ) = K/(2n£) the average superfluid 
velocity. 

Using Schwarz's metho d 35 i 36 , the three-dimensional calculation of the parti- 
cle motion in the vortex tangle was performed in the periodic box. The numer- 
ical technique was described by Samuels et al. ' . The motion, governed by 
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Fig. 13 Histogram of the angle b between u ; , and v., (vertical axis is divided by 100). From 
Kivotides, Sergeev, and Barenghi, Phys. Fluids, 20, 055105, (2008). Reprinted by permission, 
©2008 American Institute of Physics. 



Eqs. ( fl6] l. of micron-size, neutrally buoyant particle was studied in the statisti- 
cally steady state of the vortex tangle. In the three-dimensional case there is no 
need to explicitly introduce a force-free region for distances r p < a c from the vor- 
tex core: in Biot-Savart calculations, the normalization of the velocity when the 
solid particle approaches too close to a vortex is achieved by the numerical cut- 
off of the pressure gradient force acting on the particle, and a force-free region is 
automatically provided by the discretization along the vortex filament. 

Fig. 1 121 (right') shows the typical trajectory, projected on the (x, y)-plane, of 
the solid particle. This trajectory has the same features as the two-dimensional 
trajectory. Likewise, the phenomenon is observed as well of particle velocity sat- 
uration at values of (u p ) much higher than (v s ). In the three-dimensional case the 
calculated average particle velocity also agrees with scaling (124b . 

These results, together with the histogram, illustrated by Fig. [13] of the angle 
b between u p and the superfluid velocity \ s suggest that from the point of view of 
flow visualization, in the low temperature limit (at T < 0.5 K) the trajectories of 
solid particles do not reveal flow patterns of the superfluid. 



3 Self-consistent, two-way coupling model of particle-vortex interactions. 
Particle trapping on quantized vortices 

The problem of particle trapping on quantized vortices constitutes a part of a wider 
problem of reconnections of quantized vortices with the surface of the particle 
moving in the fluid velocity field. This problem cannot be analyzed based on the 
one-way coupling model introduced in Sec. [2] Instead, a more elaborate, two-way 
coupling model was developed by Kivotides, Barenghi, and Sergee v 33 1 3 V i 41 i 42 i 43 
based on dynamically self-consistent calculations which take into account an in- 
fluence of the flow field around the sphere on the evolution of the superfluid vortex. 
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3.1 Mathematical formulation 

The evolution of the vortex filament represented as a space curve X(s, f), where 5 
is the arclength, is governed by the equation 

dX/dt = \ s +\ b + \^+\ f . (25) 
The contributions are: V s - Biot-Savart integral given by formula d2 lb ; the po- 
tential field V b describes the deformation of vortices due to the presence of a 

stationary particle (on the particle surface (V s + V ) • n = 0, where n is the nor- 
mal unit vector); V 1 ^ = V</> is the potential flow field induced by the motion of the 
spherical particle: 

1 at 

<j)(x,t\r p ) = -- - _' u p -{x-r p ), (26) 

Z |X Fp| 

r 

where r p (t) is the current position of the particle centre; the contribution V is due 
to the mutual friction between the superfluid and the normal fluid: 

V f = h t , (V s + V b + ) + h*X' x (v„ - V s - V b - V ) + h„X' x(X'xv„). 

(27) 

These equations must be considered together with the equations of motion of 
neutrally buoyant spherical particle; these equations are dr p /dt = u p (f ), and 

4 3 du p <: I \ 

-na p p = 6na p ii n {Xn-Up) 

+ 2n Ps a f YS{ ^ t) + l - Ps j dS\\ s +\ h \ 2 n, (28) 

s 

where p„ is given by the first of relations dT3l) . At temperatures 1 K < T < 7\ = 
2. 17168 K, H„(T) is the viscosity of the normal fluid. At T < IK, the coefficient 
fX„ is determined by the drag force due to ballistic scattering of quasiparticles 
(phonons and rotons) off the particle surface (see the discussion later in Sec. l3.3b . 

Note that the key difference between the one-way coupling model represented 
by Eqs. dl U- (ll2b and the self-consistent, two-way coupling model considered in 
this Section is the presence of the particle-vortex interaction force represented by 
the last term in the equation J28I) of particle motion, the latter being coupled with 
equations (125l >-(127l> governing the evolution of the vortex filament. 

Numerical method of solution of the system of equations d2TT) and J25b - (l28b is 
described in detail by Kivotides, Barenghi, and Sergeev 39 (this method is based on 
generalization of the approach by Schwarz 35,44 , later developed further by Tsub- 
ota and Maekawa— , to the problem of reconnection of the quantized vortex with 
the stationary surface or the surface moving with prescribed velocity). 

Using this model we may analyze first how incorrect is the one-way coupling 
model. Fig. [14] shows the distance between the particle and the initially straight 
vortex vs time found from (a) analytical calculation 18 based on the one-way cou- 
pling model, and (b) numerical, dynamically self-consistent, two-way coupling 
model 40 given by Eqs. (121b , (125 b - d28b - As can be seen, despite the simplicity of 
the one-way coupling model, agreement is good until the moment when the parti- 
cle and the vortex are so close that the vortex reconnects to the particle surface. 
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Fig. 14 Distance between the particle and the vortex core vs time: (a) one-way coupling model; 
(b) self-consistent two-way coupling model 40 . From Barenghi, Kivotides, and Sergeev, J. Low 
Temp. Phys., 148, 293, (2007). Reprinted by permission, ©2007 Springer. 



3.2 Mechanism of particle-vortex interaction 

We will focus here only on the most important aspects of particle-vortex interac- 
tion, referring the interested reader to original publications of Kivotides, Barenghi, 
and Sergeev^ 39 . 40 ! 41 ! 42 ! 43 . 

All the calculations illustrated below were performed for the neutrally buoyant 
particle of radius a p = 1 fim located initially at the distance 2a p from the initially 
straight vortex filament. To simplify the analysis, in all examples considered below 
it was assumed that v„ = 0. 

We illustrate first an influence of the initial velocity of solid particle on the 
particle-vortex collision 42 . Fig.[T5]shows a sequence of particle-vortex configurations 
at T — 1.3K for initial velocity u„ = 25cm/s. The particle arrives from the right; 
top left: the vortex is deformed as it tries to avoid the incoming particle; top mid- 
dle and top right: the reconnection of vortex to the particle surface excites Kelvin 
waves; bottom left: the particle drags the vortex, forcing its two strands to come 
together, thus facilitating a second reconnection; bottom middle: following the 
reconnection, the vortex recoils and the particle breaks free; bottom right: the fol- 
lowing relaxation creates more Kelvin waves. 

The second computation is carried out at the same temperature for the smaller 
initial velocity u p = 20cm/s. In this case, illustrated by Fig. [16] the particle is 
trapped by the vortex: after reconnecting with the vortex the particle slows down 
and stops since it lacks the kinetic energy to stretch the vortex and to induce the 
second reconnection. Results illustrated by Figs.ll5landll6lindicate that, for each 
temperature, there exists a critical velocity, v cr of the particle-vortex approach; 
provided the relative velocity of the particle and the vortex is smaller than v cr , the 
particle will be trapped on the quantized vortex core. 

Note that in these two, somewhat artificial examples a possibility of nucle- 
ation of quantized vortices (and hence extra dissipation) by a moving particle has 
been ignored. A rather high particle velocities (20 and 25cm/s) were used for the 
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Fig. 15 (Color online) Particle-vortex collision 42 at T = 1.3K. Initial velocity u p = 25cm/s. 
Particle moves from the right. From Kivotides, Barenghi, and Sergeev, Phys. Rev. B, 77, 014527, 
(2008). Reprinted by permission, ©2008 American Physical Society. 




Fig. 16 (Color online) Particle trapping on the vortex core 4 -. T = 1.3K, initial velocity u p = 
20cm/s. From Kivotides, Barenghi, and Sergeev, Phys. Rev. B, 77, 014527, (2008). Reprinted 
by permission, ©2008 American Physical Society. 



purpose of illustration only; similar results were obtained as well for considerably 
smaller velocities. It is unlikely that in the real turbulent 4 He the particle velocity 
relative to thge vortex core can be as high as in these illustrations. For example, 
in the counterflow turbulence u p can hardly be larger than v ns , the latter usually 
being considerably smaller than the critical nucleation velocity. 

It was also found 41 that having been trapped by the quantized vortex the parti- 
cle may then drift along the vortex filament. For a neutrally buoyant, micron-size 
particle, a typical drift velocity, v drift was found to be about 0.5 cm/s. Such a drift 
can be explained by the interaction of the particle with Kelvin waves (which are 
not necessarily symmetric with respect to the particle) induced by the particle- 
vortex collision. The drift velocity provides a rather simple way of estimating the 
amplitude of collision-induced Kelvin waves, e.g. by modeling a Kelvin wave of 
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Fig. 17 (Color online) Particle trapping on the vortex core— at T = 1.3K and u p = 0. From 
Kivotides, Sergeev, and Barenghi, Phys. Fluids, 20, 055105, (2008). Reprinted by permission, 
©2008 American Institute of Physics. 

amplitude A as a vortex ring of radius A and then balancing the momentum of the 
ring with that of the drifting particle. For v drift ~ 0.5 cm/s such a procedure yields 
A 0.25 jUm. 



3.3 Influence of temperature on particle-vortex collision and particle trapping. 

Another computation- 3 - 3 - was carried out at temperature T = 1.3K for the particle 
initially at rest. The particle starts moving under the influence of the radial pressure 
gradient generated by the vortex. Fig. [T7] shows the sequence of particle-vortex 
configurations at times f = (left), 0.125 x 10~ 2 s (center), and 0.263 x 10~ 2 s 
(right). The last frame shows that the vortex traps the particle and emits a small 
vortex ring (note that this feature is not common for all trapping events) thus 
reducing the total energy of the particle-vortex configuration. The results of cal- 
culation— for the same initial configuration and particle velocity, but fi„ assumed 
to be only 0.2 of its value at T = 1 .3 K seem to suggest that at any, however small, 
non-zero initial velocity, the particle, although undergoing the process of recon- 
nection with the vortex filament, eventually breaks free, so that trapping does not 
occur. Similar (in fact, almost identical) scenario is typical of T — ► limit when 
the (viscous) damping force acting on the particle can be neglected. (However, 
this issue is less trivial than it seems - see the discussion below.) 

These examples show that the the presence of the damping is crucial for trap- 
ping of solid particles by quantized vortices, so that it can be, rather naively, ex- 
pected that trapping cannot occur at temperatures below 1 K when the normal fluid 
is absent. However, at T < 1 K there still exists a damping force, = —67ta p [l n xip 
caused by ballistic scattering of quasiparticles (phonons and rotons) off the par- 
ticle surface. Note the Stokesian form of this force (cf. Eq. (7), although n n (T) 
should now be understood not as a viscosity but as the damping coefficient). 

Calculations 33 for various values of the damping coefficient showed that trap- 
ping does not occur in the case where the damping coefficient, \i„ is smaller 
than 0.2 x ju„(1.3K). The value of the damping coefficient was measured exper- 
imentally by Jager, Schuderer, and Schoepe 46 for the spherical particle of radius 
100 jxm. Fig.ll8lshows the experimental results 46 for the coefficient X = 6na p n n 
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Fig. 18 Damping coefficient X as a function of temperature 46 (experimental results of Jager, 
Schuderer, and Schoepe, Phys. Rev. Lett., 74, 566, (1995)). Reprinted by permission, ©1995 
American Physical Society. 



in the range of temperatures including both the region T > 1 K corresponding to 
the classical viscous dissipation in the normal fluid, and T < 1 K corresponding to 
the regime of ballistic scattering of quasiparticles. It can be seen that in the tem- 
perature interval 0.6 K < T < 1 K the value of jj. n may even exceed the viscosity 
at temperatures above IK. Only at temperatures below 0.5 K the damping coef- 
ficient becomes smaller than 0.2 x jU„(1.3K). Therefore, it can be expected that 
the trapping of neutrally buoyant, 1 /xm particles does not occur at temperatures 
below 0.5 K. 

At these, very low temperatures the motion of solid particle can be modeled 
by a simpler, one-way coupling model which is based on the assumption that 
particles are not trapped on quantized vortex lines. However, at these temperatures 
the presence and motion of the particles still affects the vortex filaments, so that a 
certain modification should be necessary of the one-way coupling model. Such a 
modification was already discussed above in Sec. 12.5.3] 

It was already mentioned rather briefly that the results discussed above can be 
invalidated in the case where the extra dissipation is produced by the nucleation 
of quantized vortices in the vicinity of the particle whose velocity relative to the 
superfluid component is sufficiently large. The details of this mechanism are not 
properly understood yet and require further, rather complicated numerical study. 
Perhaps the only example of such a study so far is the work of Hanninen, Tsubota, 
and Vinen 47 who hinted at the possibility of formation of a wake, growing with 
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Fig. 19 Example 31 of the particle- vortex interaction at T = 0. Initial radius and velocity of the 
ring are R = 1.25 x 10~ 3 cm and Vr = 0.771 cm/s. The particle is initially at rest. Four frames, 
from left to right, correspond to times t = 0, 3.407 x 10~ 4 , 5.679 x 10~ 4 , and 7.698 x 10~ 4 s. 
From Kivotides and Wilkin, J. Fluid Meek, 605, pp. 367-387, (2008). ©2008 Cambridge Uni- 
versity Press. 



time, of quantized vorticity behind an oscillating sphere. However, we anticipate 
that relatively small particle velocities typical of PIV or particle tracking exper- 
iments allow to ignore the extra dissipation due to the nucleation of quantized 
vortices. 



3.4 Self-consistent model of particle collisions with vortex rings 



The self-consistent, two-way coupling model described in this Section, was re- 
cently applied by Kivotides and Wilkin 31 for numerical study of interactions be- 
tween neutrally buoyant solid particles and quantized vortex rings in the range of 
temperatures between T = and T = . It was found that trapping of particles by 
sufficiently small vortex rings never occurs, and that, at T = 0, the dominant dy- 
namical process in the particle-ring interaction is the excitation and propagation of 
Kelvin waves along the vortex ring. The collision between the particle initially at 
rest and the vortex ring of radius 1 .25 x 10~ 3 cm at temperature T = is illustrated 
by the sequence shown in Fig. [19] It was found that typical of the particle-vortex 
collisions at T = is spiraling of the particle out of the point of initial contact 
with the vortex ring. At finite temperatures the particle-vortex collision induces 
particle oscillations in the direction normal to the particle trajectory. As should 
be expected, at finite temperatures the mutual friction damps the amplitude and 
reduces the frequency of Kelvin waves propagating along the ring. 
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Fig. 20 PIV measurement of Zhang and Van Sciver 5 of v pa = u p + v s u p for particles sedimenting 
in turbulent thermal counterflow. Solid line: v pa calculated according to Eqs. 429b and < f30l > ; 
dashed line is calculated 22 from Eqs. 04t and f35t with j8 = 3. Reprinted, by permission, from 
Sergeev, Barenghi, and Kivotides, Phys. Rev. B, 74, 184506, (2006). ©2006 American Physical 
Society. 

4 Visualization experiments and their theoretical interpretation 

4. 1 Particle motion in turbulent thermal counterflow 
4.1.1 Experiment 

One of the first and, in our view, one of the most important experiments illus- 
trating strong interactions between solid particles and quantized vortices was per- 
formed by Zhang and Van Sciver 5 who studied the sedimentation of heavy (p p — 
l.lg/cm 3 ) particles in thermal counterflow produced by the heat source situated 
at the bottom of vertical apparatus, so that the normal fluid flows upwards. The ex- 
periments were performed in the temperature range from 1.62 to 2.0 K; the applied 
heat flux ranged from 1 10 to 1370mW/cm 2 . 

It seemed natural to expect that the dominating force acting on the particle will 
be the viscous drag force and, therefore, the particle velocity will be 



where v s u p is the terminal velocity of particle sedimentation given by relation 



Were this correct, the experimental data for v pa = u p + v s n„ plotted against v n 
would have collapsed on the straight solid line shown in Fig. |20] However, as can 
be seen from this Fig. [20] the results of PIV measurement showed much lower 
particle velocity. Zhang and Van Sciver 5 found that, instead of ( 130b , the particle 
velocity can be represented as u p = v„ — v s n p — v a dd, where the additional velocity, 



u p — V„ — V s ii p , 



(29) 




(30) 
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Vadd can be explained only by strong interactions between sedimenting particles 
and the vortex tangle. Zhang and Van Sciver also found that v pa Jv n ~ 0.5 inde- 
pendently of temperature, and v a dd ~ Q = pSTv n . 

4.1.2 Phenomenological theory of particle motion 

Below we will discuss the phenomenological theory, developed by Sergeev, Barenghi 
and Kivotides 22 , of the motion of micron-size particles in thermal counterflow and 
hence explain the surprising result of Zhang and Van Server's experiment. The 
relatively simple analytical model of Sergeev, Barenghi and Kivotides arises from 
the physical insight acquired using the numerical (and computationally expensive) 
two-way coupling model. 

Imagine that two strands of the quantized vortex are attached to the surface 
of spherical particle as shown in Fig. [3] (left) (at the point of reconnection the 
vortex strand is necessarily orthogonal to the particle surface). The force exerted 
by a vortex strand attached to the surface is F = J s pndS which, as shown by 
Schwarz 44 , can be written as 



The contribution of V to this force can be neglected in the case where the radius 
of curvature of the vortex strand is much larger than the particle radius. 

The vortex tangle in the counterflow can be so dense that several vortex strand 
can be simultaneously attached to the particle. Since the quantum of circulation 
is small, the leading contribution to the integral ( 131b is provided by a small area 
around the point where the vortex attaches to the surface. This enables us to find 
the following analytic approximation for the force exerted on the particle: 



It can be noticed that, in agreement with the experimental results 5 , F is a body 
force. 

Several possible particle-vortex configurations shown in Fig.|2TJcan be imag- 
ined. If configuration is symmetric, as in the first two figures, the net force is zero. 
If one or more vortex loops are asymmetrically attached to the sphere, as in the 
third and the fourth figures, the contributions from individual vortices will not 
cancel out resulting in a net body force. 

The following scenario seems realistic: the sphere, as it moves between vortex 
lines, carries along one or more vortex lines or even separate loops as the result of 
previous close encounter with vortices. This scenario corresponds to asymmetric 
particle-vortex configurations, shown in Fig.[2T] leading to emergence of the force 
exerted on the particle in the direction opposite to its motion. In the case of mod- 
erately dense tangle (a p < £) typical of the experiment of Zhang and Van Sciver 5 , 
the average force exerted on the particle by the vortex tangle can be calculated 22 
as 




(31) 



s 



F 




a N 
ln-jf£n ; . 

5 i=i 



(32) 




(33) 
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Fig. 21 (Color online) Possible particle- vortex configurations 22 . From Sergeev, Barenghi, and 
Kivotides, Phys. Rev. B, 75, 019904(E), (2007). Reprinted by permission, ©2007 American 
Physical Society. 



Here the only unknown quantity is the parameter /3 which should be determined 
by geometrical properties of the vortex tangle in the vicinity of the particle. It can 
be expected that j8 = 0(1). 

The additional velocity can now be calculated as v add = F / '(6na p jj, n ), so that 
for v pa = Up + v s up we find 

Vpa = v B - v add *[\ ^ ) v n , (34) 

where y{T) = L 1 ' 2 /v m = p s L x l 2 /(pv„) is the known mutual friction coefficient, 
with v m = v„ — v s . In agreement with the experimental results 5 we find that v pa is 
proportional to the normal fluid velocity. Moreover, the temperature dependence 
of the slope v pa /v n turns out to be the same as in the cited experiment. 

Using formula (134b . the ratio of the additional velocity to the heat flux can be 
calculated as a function of temperature: 

v add _ pK 2 ln(a p /$) f J , ( , 5) 



q \2k 2 \h„ST t 

so that Vpa can now be calculated explicitly. The dashed line in Fig. I20lreproduces 
Vpa calculated by Sergeev, Barenghi, and Kivotides 22 using formulae (I34b - d35l ) 
with j3 = 3. As can be seen, the developed phenomenological theory agrees well, 
not only qualitatively but also quantitatively, with PIV measurements of Zhang 
and Van Sciver and, therefore, seems to explain the mechanism of particle motion 
in the thermal counterflow. 

Calculations similar to those leading to formulae d33l-d35l) can also be per- 
formed in the limit of the very dense vortex tangle such that a p ^> t. In this case 
reconnections between the particle and vortices happen all the time, so that the 
particle is always attached to several vortex filaments. Calculation of the average 
force exerted on the particle by the vortex tangle yields 2 ^: 



where f5 d = 0(1) is again a geometrical factor, and (a p /£) 2 represents the cross- 
section of the particle interaction with the network of vortices. Proceeding as in 
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the previous case, a different (cf. Eq. J34b ) dependence of v pa on v„, T and a p is 
predicted: 

j3 d a p (Kyp) 2 ln(a p /%) 



" pa 



I2n 2 ix n p s f 



(37) 



where / = 1 + fyy/ (37t[i n ). Finding whether or not this prediction agrees with 
observations would require a new experiment similar to that of Zhang and Van 
Sciver 5 but for considerably higher values of the vortex line density (i.e. such that 



L > a" 2 ). 



4.1.3 Self-consistent model of particle motion in the thermal counterflow 

In order to justify the phenomenological theory of Sec. 14.1 .21 and investigate the 
particle motion in more detail, Kivotides 48,49 applied the two-way coupling, self- 
consistent model described in Sec. 13. 1 1 for a numerical study of particle interac- 
tions with the vortex tangle in thermal counterflow. Although in the cited works 
the motion of not a heavy but neutrally buoyant particle was studied, the results 
reveal some important aspects of the experiments performed by Zhang and Van 
Sciver 5 and Paoletti et al.—. 

We start this Section with reviewing the first of these papers. In calculations^ 8 - 
the normal fluid velocity in the counterflow was assumed a constant, prescribed 
value (this was the only realistic option considering the complexity of the model 
represented by Eqs. d25b -( 128l >). Numerical analysis of the particle motion was per- 
formed for the statistically steady vortex tangle modeled in the same way as in the 
earlier work 32 described in Sec. l2.5.2l 

Calculations were performed for three temperatures, 7=1.3,1.95, and 2 . 1 7 1 K. 
It was found that the following four factors strongly affect the particle motion: 
1) stratification of the vortex tangle, 2) vortex line density, 3) the average drift of 
the tangle, and 4) the intensity of Kelvin wave cascades, induced by particle-vortex 
collisions, along vortex filaments. For example, for T = 1.95K Kivotides found 
that the tangle is strongly stratified, and an average drift of the tangle is small. 
Since the vortices expand mostly in the direction normal to that of the counter- 
flow, it was also found that in the considered case the stratification does not affect 
the average properties of the particle motion. The statistically steady particle mo- 
tion is governed by the balance of the Stokes force acting in the direction of the 
normal flow, and the average particle-vortex interaction force acting in the op- 
posite direction; the latter force was found to be proportional to the vortex line 
density and the mean relative particle-tangle velocity. 

The calculations 48 confirmed the mechanism suggested by the phenomeno- 
logical theory 22 discussed in the previous Section: a formation of vortex loops 
attached to the rear part of the particle surface causes an additional force oppo- 
site to the direction of particle motion and hence reduces the slip velocity. For 
T = 1.95K and the vortex line density L = 3.284 x 10 7 cm~ 2 Kivotides' calcula- 
tion yielded (u p ) « 0.6v„, in a very good agreement with experimental results of 
Zhang and Van Sciver—. 

At the same temperature but higher vortex line density, L = 8.284 x 10 7 cm -2 
Kivotides found that the "head-on" particle-vortex collisions are more important 
and counterbalance the force caused by the formation of vortex loops in the rear 
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Fig. 22 Left: particle velocity vs time in the absence of particle-vortex collisions (the transition 
part of the curve corresponds to the particle response time). Right: u p vs time illustrating the 
particle-vortex collision; D indicates the moment when the vortex detaches from the particle 
surface and the particle velocity starts adjusting to the normal flow 49 . Adapted from Kivotides, 
Phys. Rev. B, 78, 224501, (2008). By permission, ©2008 American Physical Society. 



part of the particle, so that (u p ) fa v„. It was also found that the drift of the tangle 
in the direction opposite to the particle velocity increases the frequency of the 
"head-on" particle-vortex collisions and hence the average particle velocity. This 
effect becomes more pronounced at higher temperature; thus, for T = 2.171K 
and L = 5.846 x 10 7 cm -2 it was found that (u p ) fa 1.2v„, again in a very good 
agreement with experimental results^. 

At lower temperature, T = 1.3K, Kelvin waves generated by particle- vortex 
collisions decay slower than at higher temperatures. This, in turn, leads to large 
amplitudes of particle velocity fluctuations which become comparable with (u p ). 
The intensity of these fluctuations was found to be inversely proportional to the 
temperature and the average velocity of the tangle drift relatively to the particle. 

As argued in Sec. l3.2l (for details see original publications of Kivotides, Barenghi, 
and Sergee v 39 ' 42 ), at 0.5 K < T < Tx there should exist a critical velocity of the 
particle-vortex approach below which a neutrally buoyant particle will necessar- 
ily be trapped by the initially straight quantized vortex. Kivotides noted 49 that at 
high counterflow velocities the tangle is very dense and, therefore, particle-vortex 
interactions become so strong that it may no longer be possible to measure the 
normal velocity. In the cited work he estimated the parameters of counterflow in 
which the normal velocity would be above the trapping limit but yet sufficiently 
low, such that the tangle is sufficiently dilute to enable the measurement of the 
normal velocity by the PIV or the particle tracking technique. 

Using the same model as in his previous work 48 , Kivotides 49 analyzed nu- 
merically, at temperature T = 1.3K, the particle motion in the thermal coun- 
terflow with the normal velocity v„ = lOcm/s and the vortex line density L = 
1.168 x 10 6 cm~ 2 . Performing calculations for different initial positions of the 
particle, Kivotides found that in more than 50% of realizations the particle moved 
through the computational domain without experiencing any collision with the 
vortex tangle. On average, the deviation, caused by particle-vortex collisions, of 
the particle velocity from that of the normal fluid was found to be less than 4%. 
Moreover, in the case where the particle does not follow the normal fluid it tracks 
the motion of the vortex tangle. On the time series of particle velocity the events of 
the particle-vortex collision can be identified by strong oscillations (Fig. |22] (right), 
cf. the left frame showing the velocity of the particle moving through the tangle 
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Fig. 23 (Color online) PDFs of the particle velocity component in the direction of counter- 
flow. After Paoletti, Fiorito, Sreenivasan, and Lathrop, /. Phys. Soc. Japan, 77, 1 1 1007, (2008). 
Reprinted by permission, ©2008 Physical Society of Japan. 

without collisions with quantized vortices) which can be filtered out to restore the 
normal velocity. 

4.1.4 Particle tracking experiments 

Another important visualization experiments in thermal counterflow are those by 
Paoletti, Fiorito, Sreenivasan, and Lathrop 10 . In contrast with the PIV technique, 
which analyzes the local average properties of the particulate flow, the parti- 
cle tracking technique investigates individual particle trajectories. Experiments^ 
were performed in the vertical apparatus with the heater at the bottom. The tem- 
perature and the applied heat flux, q ranged from 1.8 to 2.15K and from 13 to 
91mW/cm 2 , respectively (cf. the experiment of Zhang and Van Sciver- with q 
between 110 and 1370mW/cm 2 ). Tracers were solid, micron-size hydrogen par- 
ticles of density slightly smaller than that of liquid helium. 

It was found that two distinct types of particle trajectories can be identified: 
1) smooth trajectories corresponding to particles moving upward in the direction 
of the normal flow, and 2) irregular trajectories of particles moving downward; 
the latter trajectories were those of particles trapped on quantized vortices. Exper- 
imental observations 10 seem to confirm the earlier theoretical results™ (see above 
Sec. 13.2b that at sufficiently low relative velocities between the normal fluid and 
quantized vortices the particles should be trapped more easily by the vortex tan- 
gle; at high relative velocities, although the particles interact with the tangle, it is 
unlikely that they become permanently trapped on the vortex lines. 

Calculated from the experimental data 10 , the probability distribution function 
(PDF) of the velocity component in the direction of the counterflow is bimodal, see 
Fig.|23j with the right peak at the normal fluid velocity, and the left, much broader 
peak corresponding to particles trapped by the vortex tangle. Note that, as would 
be expected, the fraction of particles trapped by the tangle and hence bimodality of 
the PDF disappear either with increasing temperature at constant heat flux (so that 
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the vortex line density becomes smaller), or with increasing heat flux at constant 
temperature (so that v„ becomes sufficiently high to prevent particle trapping on 
quantized vortices). 

Discussing these experimental findings, Paoletti et al. 10 claimed that their re- 
sults do not agree with observations of Zhang and Van Sciver 5 who found that 
the particle velocity is proportional to the normal velocity, u p ps 0.5v„ indepen- 
dently of temperature. The authors 10 attributed this discrepancy to Zhang and Van 
Sciver's PIV technique which measures the local average properties of the partic- 
ulate flow, while the technique employed by Paoletti et al. tracks the motion of 
individual particles. Paoletti et al. also stressed that their interpretation of exper- 
imental results differs significantly from the theoretical explanation of Sergeev, 
Barenghi, and Kivotides 22 (see also Sec. 14.1 .21 above) whose underlying assump- 
tion was that every particle is affected by quantized vortices as it moves through 
the tangle, while the experimental observations 10 showed that there is a signif- 
icant fraction of particles which move freely through the tangle without experi- 
encing particle-vortex collisions. Furthermore, Paoletti et al. observed a signif- 
icant temperature dependence of the particle motion, while the experiment 5 of 
Zhang and Van Sciver and the phenomenological theory 22 of Sergeev, Barenghi, 
and Kivotides both show that the disparity between the particle velocity and v„ is 
practically independent of temperature. 

Below we will argue that the contradiction between the experimental results 
of Zhang and Van Sciver 5 and Paoletti et al.—, as well as between theoretical 
interpretations of these results by Sergeev et al. and Paoletti et al. is only apparent. 
In fact, it seems that two recently published works by Kivotides 4 -^^, reviewed 
in Sec. 14.1.31 already resolve this apparent contradiction. In these works, based 
on the two-way coupling, self-consistent (practically "first-principle") approach, 
Kivotides has shown that there exist two distinct regimes of particle motion, one 
corresponding to a relatively dense, and another to a relatively dilute vortex tangle. 
In the first of these publications Kivotides found that in the case where the vortex 
tangle is sufficiently dense (L = 3.284 x 10 7 cm -2 corresponding to the intervortex 
spacing t w 1.7/Jm in the considered example 48 ) the average particle velocity 
is 0.6v„ in a very good agreement with the experimental results of Zhang and 
Van Sciver 5 . Kivotides also showed that in the considered case the particle is 
permanently affected by quantized vortices as it moves through the tangle, and 
that the mechanism of particle-vortex interactions agrees with that suggested by 
the phenomenological model of Sergeev, Barenghi, and Kivotides 22 . On the other 
hand, in the second of his publications 49 Kivotides, having analyzed numerically 
a particle motion in a more dilute tangle (L = 1 . 168 x 10 6 cm -2 , 1 sa 9 jUm), found 
that in more than half realizations the particle moves through the tangle with the 
normal fluid without interacting with vortices. In the remaining less than 50% 
realizations he observed strong particle-vortex interactions which in most cases 
can be described as trapping-untrapping events. Were it calculated based on the 
results reported in the second of his publications 49 , the PDF would have the same 
bimodal shape as found experimentally by Paoletti et al.—. 

Although, in terms of the heat flux, the regimes of counterflow in two reviewed 
experiments were adjacent, in most observations of Zhang and Van Sciver the 
heat flux was an order of magnitude or more higher than that in the experiments 
of Paoletti et al. This means that, at the same temperature (e.g. 1.95K in both 
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experiments), the intervortex spacing in the experiments of Zhang and Van Sciver 
(I ~ 6/im with 1.7/im particle in a typical experiment) was at least an order of 
magnitude smaller than in the reviewed experiment of Paoletti et al. Therefore 
it can be argued that two reviewed experimental observations do not contradict 
each other but simply correspond to two distinct regimes of particle motion. The 
parameter defining each of these regimes is the ratio of the particle size to the 
intervortex spacing, a p /t. 

As far as the issue of temperature independence of the factor k in the relation 
(u p ) ~ kv n is concerned, Zhang and Van Sciver's results 5 are not truly temper- 
ature independent: a closer inspection reveals a relatively weak dependence of k 
on temperature (the value k sa 0.5 was obtained by averaging of a large number 
of experimental data). Likewise, in the phenomenological theory developed by 
Sergeev, Barenghi, and Kivotides 22 this factor is only approximately temperature 
independent being a function of the mutual friction coefficient y and of the param- 
eter P which, characterizing the geometry of interactions between the tangle and 
the particle, is itself temperature-dependent. 

4.1.5 P IV experiment in thermal counter/low with cylindrical obstacle 

Among surprising experimental results obtained by the PIV technique is the re- 
cent observation by Zhang and Van Sciver 6 of the apparently stationary normal 
fluid eddies in the thermal counterflow past a cylinder. In the cited work, Zhang 
and Van Sciver visualized the motion of small particles in the thermal counterflow 
around the cylinder of diameter D = 0.635 cm fixed in the center of rectangu- 
lar channel of a cross-section 3.89 x 1.95cm 2 . The counterflow was produced, in 
two separate experiments, by the heat flux q = 0.4 and 1.12W/cm 2 at tempera- 
tures T = 1.6 and 2.03 K, respectively (corresponding to the Reynolds numbers 
Re = pDv n /n n = 4.1 x 10 4 and 2.1 x 10 4 ). Solid particles used for visualiza- 
tion in the PIV experiments were polymer microspheres of diameter 1.7/im and 
density 1.1 g/cm 3 . In these experiments Zhang and Van Sciver observed the for- 
mation of large-scale eddies of the particulate motion located both downstream 
and, surprisingly, upstream of the cylinder with respect to the normal flow. These, 
apparently stable vortices of the particulate flow field were located at distances 
about 3 cylinder radii from its center at the angles ±45° and ±135° to the axis 
along the undisturbed flow through the center of the cylinder, see Fig. [24] Note 
that the observed flow structures do not have a classical analogue. 

Zhang and Van Sciver attributed the existence of apparently stationary nor- 
mal eddies to the mutual friction interaction between quantized vortices and the 
normal fluid. However, our recent study 50 showed that perhaps the experimental 
results 6 can be interpreted without invoking the mechanism of interaction between 
the normal fluid and quantized vortices. Indeed, the calculation of motion of point 
vortices in the imposed potential flow around the circular disk shows that there 
exist stationary locations of point vortices, both at the rear and at the front of the 
disk. These locations are unstable: any perturbation of the initial stationary posi- 
tions of point vortices leads, eventually, to sweeping of point vortices away from 
their initial locations. Furthermore, some of these stationary locations are posi- 
tioned practically as the eddies seen by Zhang and Van Sciver. The point vortices 
in the vicinity of such positions will remain close to their initial locations during 
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Fig. 24 (Color online) Streamlines of the particulate motion in the counterflow around the cylin- 
der^ (q = 1 . 12W/cm 2 , T = 2.03 K). From Zhang and Van Sciver, Nature Physics, 1, 36, (2005). 
Reprinted by permission, ©2005 Macmillan Publishers Ltd. 



the time period corresponding to the duration of the experiment 6 and hence seen 
as apparently stable. 

Although this, purely classical explanation of the apparent stability of vortex 
structures does not invoke any interaction between the normal and the superfluid 
components of 4 He, the emergence of eddies seen in the experiment 6 might still 
require an explanation based on the analysis of mutual friction between quantized 
vortices and the normal fluid. 



5 4 He channel flow and turbulent boundary layer 

This short Section describes the recent experiment which might be a beginning 
of systematic study of nonuniform 4 He flows. Xu and Van Sciver 51 reported the 
PIV measurements, using micron-size deuterium particles, of the 4 He forced flow 
in the rectangular channel for the normal fluid Reynolds numbers ranging from 
9 x 10 4 to 4.5 x 10 5 and temperatures from 1.65 to 2.10K. In this experiment, at 
scales larger than the intervortex spacing, the normal and superfluid components 
of 4 He can be considered as fully interlocked so that the measurements of the 
particulate velocity field yield an unambiguous velocity profile of the fluid. The 
results of Xu and Van Sciver's measurements were summarized in two graphs^- 
shown in Fig. 1251 

Xu and Van Sciver concluded 5 1 that in the wall region the velocity distribution 
agrees reasonably well with the classical nth-power law (for n ranging from 7 to 
8.8). They also addressed the following questions which yet to be answered. Why 
the velocity profile is wider and flatter than that in the classical viscous channel 
flow? What is the nature of dependence on Reynolds number of the results shown 
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Fig. 25 (Color online) Normalized velocity profile 51 for different Reynolds numbers at T = 
2.10K (left), and at different temperatures for the mean velocity U mea „* = 11.5cm/s; W is the 
width of the square channel. From Xu and Van Sciver, Physics of Fluids, 19, 071703, (2007). 
Reprinted by permission, ©2007 American Institute of Physics. 

in Fig. [25] (left)? Why, as seen from Fig. [25] (right), the normal fluid density does 
not seem to affect the shape of velocity profiles? 

6 Visualization of vortex reconnections. Velocity statistics in decaying 
quantum turbulence 

In a recent experiment 9 Bewley, Paoletti, Sreenivasan, and Lathrop observed the 
motion of solid hydrogen particles trapped on quantized vortices. The specific 
purpose of this work was a direct experimental investigation of vortex reconnec- 
tions in turbulent 4 He. A sequence of images 9 illustrating the motion of particles 
trapped on vortex filaments is reproduced in Fig. 1261 

To quantify their results, Bewley et al. 9 assumed that the evolution of recon- 
necting vortices can be characterized by a single scale parameter, l(t). Using, as a 
measure of /(f), the experimentally observed distance between two particles clos- 
est to the point of reconnection, they found that the evolution obeys the scaling 
I ~ (t — to) 1 / 2 , where to corresponds to the moment of reconnection of two vor- 
tices. 

The particle tracking technique was further developed by Paoletti, Fisher, Sreeni- 
vasan, and Lathrop 11 to investigate, by analyzing the trajectories of tracer parti- 
cles, the velocity statistics in decaying quantum turbulence. The decay of turbu- 
lence produced initially by the thermal counterflow was studied after the coun- 
terflow has been stopped by switching the heater off. Paoletti et al. found that the 
PDF of the particle velocity is strongly non-Gaussian with a pronounced tail obey- 
ing v~ 3 power law, see Fig.|27j The authors attributed such a tail to high velocities 
produced by reconnections of quantized vortices. The experiment has stimulated 
calculations 52 of turbulent velocity statistics in the context of three-dimensional 
and two-dimensional atomic Bose-Einstein condensates; after finding similar non- 
Gaussian velocity statistics in these related quantum systems the authors argued 
that, in general, non-Gaussian statistics arise from the singular nature of the ve- 
locity field around quantized vortices. 

The results of fundamental importance reported in these two reviewed publi- 
cations^^ raise, however, the following question: what effect have the particles 
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Fig. 26 (Color online) After Bewley, Paoletti, Sreenivasan, and Lathrop (Proc. Nat. Acad. Sci., 
105, 13707, (2008)): visualization 9 of the vortex dynamics. The particles trapped on quantized 
vortices can be easily identified. In each of the sequences, a, b, and c, the images were taken at 
50ms intervals, (a): two approaching vortices have several particles trapped on their filaments 
(the first frame shows the projection in which the vortices appear crossed). The sequence (b) 
seems to show the vortices moving apart after the reconnection. The sequence (c) illustrate the 
authors' 9 method of identification of reconnecting vortices by a sudden motion of two tracer 
particles away from each other. ©2008 USA National Academy of Sciences. 



trapped on vortex filaments on the evolution of quantized vortices and, in partic- 
ular, on the reconnection of two approaching vortices? In their work— reviewed 
earlier in Sec. 14.1 .41 the authors claimed that, because in their experiments the dis- 
tance between micron-size particles trapped on the same vortex filament was typ- 
ically about 100 fim, the influence of trapped particles on the evolution of vortices 
can be neglected. This claim is yet to be justified (especially where the motion 
and evolution of vortex filaments near the reconnection point is concerned). It is 
unlikely that at present this problem can be treated experimentally; a more realis- 
tic alternative seems to be a development of the mathematical model of a vortex 
filament loaded by trapped solid particles^. 
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Fig. 27 (Color online) Probability distribution functions ' 1 of the velocity components in the 
direction of (v-) and normal to (v r ) the counterflow. The distributions are scaled with o~,,_ = 
0.074 cm/s and a Vx = 0.066 cm/s, respectively. The dashed (blue) line - velocity PDF, scaled 
with rj,, = 0.25 cm/s, for classical turbulence. From Paoletti, Fisher, Sreenivasan, and Lathrop, 
Phys. Rev. Lett., 101, 154501, (2008). Reprinted by permission, ©2008 American Physical 
Society. 

7 Related techniques: Detection of vortices by trapping of negative ions and 
imaging of He2 molecules 

7.1 Ion trapping 

The technique of vortex detection by ion trapping was recently reviewed in a num- 
ber of regular research publications including some in this journal 54,55 , so that we 
will present here a rather brief overview of theoretical and experimental aspects 
of this technique. 

The idea of using charge carriers for detection of quantized vortices in super- 
fluids is much older than that of using solid particles for the flow visualization in 
4 He and dates back to the works of Careri 56 , and Reif and Meyer 57 . The technique 
is based on the phenomenon that a negative charge (a single electron) injected into 
4 He self-localizes itself in a spherical void (known as the ion or electron bubble) of 
radius 12 — 20 A (depending on pressure) from which helium atoms are excluded. 
(A relatively large size of this, almost macroscopic bubble justifies inclusion of 
this Chapter into our review otherwise concerned with the motion of much larger, 
solid particles in 4 He.) The model of the negative ion bubble was proposed by Fer- 
rel and later elaborated by Kuper and other authors 58 and confirmed experimen- 
tally by Levine and Sanders 59 . Consistent with the bubble model of the negative 
charge carrier is the effective mass of the ion, M* determined in Ref. 60 . It appears 
in the Langevin equation of motion for the ion's drift velocity, Uq in the electric 
field, E: eE = M* (dup jdt + Up /t), where e is the electron charge, and T the phe- 
nomenological relaxation time. The effective mass of the ion was found 60 to be 
of the order 100 masses of 4 He atom (in agreement with the bubble model), and 
% appears to be the viscous relaxation time introduced above in Sec. 12.21 (The de- 
tailed treatment for the Stokes problem of the ion bubble motion, albeit in normal 
4 He, was given by Ostermeier and Schwarz- 6 -^.) 
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The possibility of using negative ions as detectors (or "probe particles") of 
quantized vortices is due to the fact that the ion bubble is attracted to the vortex by 
the Bernoulli force. The combined action of the Bernoulli effect and the reduced 
condensation energy of the core produces a potential well of the depth of the 
order 50K 62 and hence traps the ion. This mechanism 62 is practically identical to 
that described above in Sec. l2.3l of trapping of solid particles on the vortex cores. 
(Note that the mechanism of trapping of the ion bubble by the quantized vortex in 
4 He modeled as the Bose condensate was studied by Berloff and Roberts 23 and 
described above in Sec. l2.3.2"l ) 

The model of the positive charge carrier in 4 He was developed by Atkins 63 
in 1959 and became since then commonly accepted. According to this model the 
positive ion exerts an electrostrictive attraction on the surrounding fluid thus caus- 
ing a liquid-solid transition resulting in a core ("cluster" or "snowball") of solid 
helium. The radius of this snowball is from 7 to 9 A, depending on pressure. Like 
the negative ion bubble, such a snowball will be attracted, by the Bernoulli force, 
to the vortex core so that, in principle, positive ions can also be used for detection 
of quantized vortices. However, a smaller than that of the ion bubble radius of the 
snowball and, therefore, the cross-section of the snowball-vortex interaction make 
the experimental technique using positive ions more difficult and less practical. 
Hence in this review we will discuss only works concerned with negative charge 
carriers. 

The realization of the negative ion technique can be illustrated by the follow- 
ing arrangements of the experiment, typical of the early observations of single 
vortex lines in rotating 4 He: quantized vortex lines are charged first by ion bub- 
bles trapped on them using e.g. the electron beam emitted orthogonally to the axis 
of rotation. The vortex lines are then detected by applying the electric field parallel 
to the axis of rotation so that trapped electron bubbles slide along vortex filaments 
to a collector attached to an electrometer. The latter registers the amount of col- 
lected charge which is proportional to the number of vortex lines present. In its 
early, simplest version this technique allowed to determine only the total number 
of (almost) straight quantized vortex lines in a slowly rotating container. Its sub- 
sequent modifications made possible also detection of individual vortex lines and 
later the measurement of the vortex line density in turbulent 4 He. 

The theory of ion trapping by the quantized vortex has been developed in clas- 
sical works of Donnelly, Roberts, and Park s 21 ' . To analyze ion-vortex interac- 
tions, two competing mechanisms were taken into account: trapping of the ion 
bubble, moving in the electric field, in the potential well of the vortex, and escape 
of the trapped ion due to its Brownian motion, which was considered to be in equi- 
librium with thermal (quasiparticle) excitations in 4 He. Two parts of this problem 
were calculations of 1) the trapping (or capture) cross-section, a, cm, and 2) the 
escape probability, P, s _1 ; the second part was analyzed using Smoluchowski or 
Fokker-Planck equations in the framework of the Kramers-Chandrasekhar method 
for calculation the probability of escape of a particle by diffusion from a poten- 
tial well. The capture cross-section, which was found to be of the order 10~ 6 
to 10 _5 cm and to decrease with the electric field and increase with temperature 
to about 1.6K, at which temperature it drops sharply. The escape probability (or, 
equivalently, the mean trapping lifetime, tg_ ~ 7 3-1 ) was found in a good agreement 
with experimental data obtained by Springett, Tanner, and Donnell y 65 ' as well 
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as with experimental results of Douglass and Cade . The effective cross-section 
can then be calculated as a e ff = ae~ Pl , where f is some characteristic time. Parks 
and Donnelly 21 further developed the theory in order to make possible calculation 
of the ion bubble radius from the experimentally measured mean trapping lifetime. 
This required, in particular, a calculation of the so-called substitution energy, i.e. 
a kinetic energy of rotating superfluid excluded by the trapped bubble. For the ion 
bubble of radius R trapped symmetrically on the vortex core this was found as 



(note that formula (fT8T l of Sec. l2.3.l1 follows from Eq. d38l l assuming R = a p S> £). 
In the context of further experimental studies discussed below, it is worth noticing 
that the substitution energy (equal to the depth of the potential well) is about 50 K. 
Parks and Donnelly's study of bubble radii was further developed by Springett and 
Donnelly 68 who used the measurements of trapping cross-sections in the rotating 
container to deduce that the radius of the ion bubble decreases with pressure. The 
theory developed by Donnelly et al. also predicted that the mean trapping life- 
time should increase with the superfluid density and with the radius of ion bubble. 
These predictions were soon confirmed experimentally by Springett 65 who mea- 
sured, at various pressures and temperatures, the cross-section of the ion capture 
and ion mobilities in order to derive from these data the pressure and temperature 
dependence of the ion bubble radius and the mean trapping lifetime. Later Pratt 
and Zimmermann 69 also measured the mean trapping lifetime in the wide range 
of temperatures and pressures (from vaporization and solidification) and showed 
that at constant pressure the lifetime rapidly decreases with T. In particular, it was 
found that at saturated vapor pressure trapping becomes negligible at tempera- 
ture above 1.7K, the value which is now often referred to as the "abrupt lifetime 
edge". It was also shown that the temperature, below which trapping becomes sig- 
nificant increases with pressure. Glaberson 70 measured the mobility, as a function 
of temperature and pressure, of negative ions trapped on quantized vortex lines, 
and arrived at the important conclusion that the negative ion bubbles do not de- 
form as they are trapped on vortex cores. He also developed a new model, which 
appeared to account satisfactory for available experimental results, for the drag 
exerted on the trapped ion bubble, based on the assumption that the vortex line 
has several Athick central core surrounded by a tail of excess roton density with 
momenta opposite to the direction of circulation. 

Among the first applications of the ion trapping technique for an investigation 
of vortex structures in 4 He was an observation by Northby and Donnelly— of a 
nonlinear dependence of the number of quantized vortex lines on the angular ve- 
locity of the rotating container. This nonlinearity was found to be associated with 
the existence of a near-wall region, observed in the cited work, free of quantized 
vortices. 

The early experiment s 56 ' 72 have indicated that ions interact strongly with tur- 
bulence in 4 He, but no quantitative data have been obtained yet. The first experi- 
mental study of the turbulent vortex tangle by the ion trapping technique was un- 
dertaken by Sitton and Moss 73 in 1969 (in fact it was the first direct experimental 
confirmation that the turbulence in the superfluid component of 4 He is composed 
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of individual, quantized vortex lines indistinguishable from those produced in the 
(slowly) rotating container except for their configurations). The vortex tangle was 
produced by the supercritical heat current at temperature above 1.6K. Since most 
of the vortex lines are no longer straight but have a configuration of loops and 
kinks, the charge can no longer slide along the filaments to a collector but is 
trapped inside the tangle. The fraction /q of the trapped charge was measured 
and then linked with the vortex line density by the relation /q = (1 +P/ui<yL)~ 1 , 
where H; is the ion velocity, and L can be linked with the heat flux and properties 
of 4 He by the well known Vinen equation—. 

A series of experimental studies and their interpretation, beginning 1972 and 
spanning the period of nearly 30 years, was undertaken by the group led by Packard 
and Williams. In 1972, Packard and Sanders 75 showed the possibility of detect- 
ing individual vortex lines in rotating 4 He. The experiments were performed in a 
cylindrical container whose rotation slowly accelerated so that the number of vor- 
tex lines present increased with time. The experimental arrangements were similar 
to those described above,. i.e. the charge trapped on vortex lines was measured by 
applying the axial electric field thus transferring the charge to a collector attached 
to an electrometer. The appearance of each new vortex line was detected by a step- 
like increase in the electrometer's reading. Bringing the container slowly to rest, 
the authors also discovered the existence of remanent vorticity studied later by 
Awschalom and Schwarz (see below). 

Later Williams and Packard 76 developed the first photographic technique which 
directly visualized spatial positions of individual quantized vortex lines in rotating 
4 He. The experimental arrangements were similar to those of Ref. 75 , but in addi- 
tion the magnetic focusing was used to stabilize beams of electrons emitted from 
the charged vortex lines; these beams impinged on a phosphor screen so that the 
positions of quantized vortices could be actually photographed. Experiments were 
performed at temperatures lower than 0.3 K, and up to 0.8% of 3 He (which acted 
as a fixed amount of normal fluid) was added to stabilize the vortices. Among other 
results, it was found, in agreement with the classical Feynman's prediction, that 
the average intervortex spacing is ~ ^/h/lconu, where (0 is the angular speed of 
rotation, but in this work an expected stable triangular lattice of vortices was never 
observed. The photographs of regular, symmetric arrays of vortices were obtained 
when this experiment was repeated some years later 77 . (These photographs should 
be very familiar to a reader through numerous reproductions in other publications 
and hence are not shown here.) The authors ascribed to mechanical disturbances 
their earlier failure to observe the regular structures. They also argued that the 
symmetric state of the system is determined not by the absolute minimum of the 
free energy, as has been believed earlier, but by the combination of the past his- 
tory of the system and the local minimum of the free energy, concluding that the 
symmetric state is highly metastable. The detailed description of the photographic 
technique for visualizing the positions of quantized vortices in rotating 4 He was 
given in work 78 . Much later one of the authors claimed 79 that usual helium ions 
are not really suitable for detection of thermally excited vortices, and that multi- 
electron bubbles should be employed for this purpose. Also such an experiment 
was outlined in the cited paper, the authors of this review are not aware of any 
practical development in this direction. 
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Using the experimental technique of Ref. , this research group also revis- 
ited 80 the problem, analyzed earlier theoretically by Donnelly, Roberts, and Parks 2 
of the lifetime of ions trapped on the vortex lines. Theoretical analysis predicted 
that in the temperature interval 0.6 < T < 1 . 1 K the mean lifetime should be longer 
than 10 13 s and that it should increase with decreasing temperature. However, the 
experiments revealed that the observed lifetime was many orders of magnitude 
smaller and was actually decreasing as temperature decreased. To resolve the con- 
tradiction with the theoretical predictions, the authors suggested that at tempera- 
tures lower than 1.5 K the trapping lifetime is no longer determined by the intrin- 
sic properties of bubble-vortex interaction but by a time scale of the vortex motion 
until it is destroyed by the container's wall. The details of this mechanism were 
suggested by the authors as follows: as the moving vortex filament encounters the 
wall, it is destroyed (so that the charge is collected by the wall), but, to maintain 
the equilibrium value of the vortex line density, another, uncharged filament is 
created. In the case where the time scale of migration of the vortex filament to 
the wall is smaller than the intrinsic trapping lifetime of the ion bubble, the vortex 
motion will be a limiting factor of the charge loss. The authors concluded that 
at low temperatures such that the normal fluid density and hence viscous damp- 
ing (by the mutual friction) become negligible, the observed lifetime is actually 
the measure of the timescale of the vortex migration to the wall and, therefore, 
should be independent of temperature, in agreement with the authors' observa- 
tions. Furthermore, the authors argued that their conclusions can also be related 
to the earlier observations by Cheng, Cromar, and Donnelly 81 that the ion trap- 
ping reduces in the presence of the axial heat current. The explanation suggested 
in Ref. 80 is that in the counterflow turbulence the mutual friction increases the 
intensity of the vortex motion and hence the rate of destruction of vortices by the 
wall. 

Employing the ion trapping technique, one of the works of key importance 
for understanding the mechanism of the onset of quantum turbulence is that of 
Awschalom and Schwarz 82 . This experiment used two parallel plates, immersed 
in 4 He, one of which serves as a charge collector detecting the vortices pinned to 
both plates. (The experiment aimed, in particular, to support the earlier Schwarz's 
idea 83 , based on his analysys of vortex reconnections, that under some condi- 
tions vortex singularities can multiply.) The experiment showed that quiescent 
4 He contains a rather large number of quantized vortex filaments which are pinned 
metastably to the parallel plates. The line density of these remanent vortices seems 
to be history independent, and its existence was established upon going, in temper- 
ature, down through the X transition. (Note that in this work the vortex line den- 
sity was not measured directly but estimated, considering the rather complicated 
geometry of electric field lines, from the measurement of the collected charge.) 
Awschalom and Schwarz found that the line density of remanent vortices is very 
close to the critical line density, introduced by Tough 84 , below which vortices 
disappear. The authors' interpretation of the critical line density was that the vor- 
tex filament do not disappear but become immobilized by pinning to the walls 
and hence cannot be observed by any conventional experimental technique. The 
authors also justified their ideas by analyzing the evolution of the vortex line den- 
sity in decaying turbulence produced by ultrasound. They observed that, after the 
ultrasound has been switched off, L decays not to zero but to the value of the rema- 
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nent vortex line density. (Earlier experiments of Milliken Schwarz, and Smith- 
produced qualitatively similar results.) The authors concluded that any volume 
"will be penetrated ab initio by quantized vortices stabilized by surface pinning", 
and hence all experiments on the evolution of the vortex line density should be 
interpreted as if L starts from the value corresponding to the remanent vortex line 
density. 

The first attempt of systematic experimental study, by means of the ion trap- 
ping technique, of the decay of quantum turbulence was made by Davis, Hendry, 
and McClintock—. In turbulence, generated by the oscillating grid, some ions 
emitted from the tip get trapped on vortices within the tangle and hence reduce 
the current arriving at the collector. The evolution of the vortex line density can be 
estimated from the evolution of this current, and the experiment clearly demon- 
strated the production and decay of quantum turbulence as well as its spatial dis- 
tribution. However, this experiment had certain disadvantages. In particular, the 
value of the vortex line density in the low temperature limit (from 200 down to 
22 mK) remained unknown, and the data on trapping cross-sections at these tem- 
perature had not yet been available. Nevertheless, thies experiment yielded a very 
useful information on the time scale of decay, and also showed that in the temper- 
ature range from 22 to 70 mK the process becomes temperature independent. 

Clearly, a modification of the classical ion trapping experimental arrangements 
was required that would make possible, in the low temperature limit T < 200 mK, 
to measure directly the dynamics of the vortex tangle as well as the trapping cross- 
sections. Such a new technique was recently developed by Walmsley, Golov and 
their co-worker s 54 i 55 - 87 i 88 ; the detailed description of the experimental cell and 
methods can be found in the first two of the cited works. 

The experimental cell, mounted on a rotating cryostat, is a cube whose sides 
are electrodes serving as the charge collectors. In the earlier experiment 54 just 
one ion emission tip was fitted at the bottom side of the cell, while later exper- 
iment 55 - 8788 ' 89 used two emission tips fitted at the bottom and one of the side 
plates. In order to provide the means of measuring the spatial properties of the 
vortex tangle, a difference of electric potential was kept between some of the 
electrodes to ensure depletion of the ion current emitted by the tips. An analy- 
sis of ion trapping by quantized vortices suggests the exponential decay of the ion 
current, collected by the electrodes, with the vortex line density, so that the lat- 
ter can be recovered from the measurements of the current I(t) from the relation 
L(t)/La = (ad)- 1 ln[/(°°) //(*)], where d is the size of the experimental cell. The 
turbulence was produced by spin up or/and spin down of the rotating experimental 
cell. 

In Ref . 54 the authors argued in favor of using, as probe particles, the charged 
vortex rings rather than bare ions. The reason is that the trapping cross-section 
of the bare ion bubble is small, a ~ 10~ 6 cm 64 - 90 , and was never measured at 
temperature below 0.8K. On the other hand, at low temperatures, in the case 
where the vortex nucleation velocity, v c is lower than Landau critical velocity, 
Vl ~ 5 x 10 3 cm/s, the ion, upon reaching the velocity v c will nucleate the vortex 
ring that captures the ion thus forming the stable ion-ring complex. The details 
of this mechanism can be found in the classical paper by Rayfield and Reif— 
(see also the work by Nancolas, Bowley, and McClintock 92 and Donnely's mono- 
graphs 62 ). While attempting, initially, to repeat the experimental conditions of the 



42 



work by Davis, Hendry, and McClintock , who tried to eliminate nucleation of 
vortex rings by ions, Walmsley, Golov and their co-authors eventually came to a 
conclusion that at low temperatures charged rings are preferable as proble parti- 
cles. The trapping cross-section of the charged ring is of the order of its diameter, 
a ~ 1 /im 93 , which is much larger than the trapping cross-section of the bare ion 
bubble, and hence it is considerably easier to detect quantized vortices by charged 
rings (as was stated much earlier by Guenin and Hess 94 ). Some disadvantage of 
using charged vortex rings is their dynamics which is more complicated than that 
of bare ion bubbles (we refer the reader to the theoretical studies^). 

The experimental studie s 54 ' 55 of Walmsley et al. showed that using charged 
rings rather than bare ions it is possible to detect the change in the vortex line den- 
sity upon starting and stopping rotation at temperatures from 0.5 down to 30 mK. 
Monitoring in Ref. 55 the trapping of charged rings by rectilinear vortices gener- 
ated by slow rotation the authors were able to measure the ring-vortex trapping 
cross-section at these temperatures. Furthermore, at low temperatures injection of 
charged rings can be used to create a vortex tangle^. 

An application of this technique to the study of quantum turbulence in the limit 
T — > has already brought some non-trivial results. Thus, Walmsley et al.— stud- 
ied, at temperatures from 1.6 down to 0.08 K, a decay of the homogeneous vortex 
tangle produced by a sudden spin down of the rotating cell. They found that the 
decay of the energy flux, £ is the same as that of the classical turbulence at high 
Reynolds numbers, i.e. £ ~ f -3 , and that the vortex line density decays as f~ 3 / 2 . 
Most importantly, it was found that at T « 0.8 K the effective kinematic viscos- 
ity, Vk drops sharply from the value Vk ~ 0.1 JC and approaches Vk ~ 0.003 k: as 
T — > 0. The authors linked this drop to the transition to the new form of turbu- 
lence, the so-called "ultraquantum" (or "Vinen" as opposed to the classical, "Kol- 
mogorov") turbulence. This, less structured than (quasi)classical regime of turbu- 
lence is characterized by the cascade in which the energy is transferred to smaller 
scales by Kelvin waves on individual vortex lines (and, therefore, by the absence 
of any large scale motion). The mechanism of energy dissipation is though to be 
acoustic: smallest perturbations are emitted as phonons. Clearly, the ultraquantum 
regime should be characterized by an effective kinematic viscosity, Vy different 
from the quasiclassical Vk- 

It should be emphasized, though, that the regime of turbulence in the cited 
work 88 was essentially quasiclassical due to forcing at large scales (rotation). The 
existence and properties of ultraquantum regime at temperatures below 0.5 K were 
addressed a year later in the recent work of Walmsley and Golov 89 . The vortex 
tangle with was produced not by spin down of the rotating scale but by charged 
vortex rings so that no large-scale flow was generated. It was found that in the 
ultraquantum regime the vortex line density decays as f _1 , and that the effective 
kinematic viscosity of the Vinen turbulence is Vy ~0.1 K. 

In conclusion of this Section we would like to emphasize that the considered 
so far two experimental techniques of visualization (or, at least, detection) of flow 
properties in turbulent 4 He, i.e PIV/particle tracking and the ion trapping tech- 
nique do not compete but rather complement each other. Indeed, the ion trapping 
technique detects quantized vortices without disturbing the flow, and it is efficient 
in the low temperature regime; at temperatures above 1.7K ion bubbles can no 
longer be trapped on quantized vortices. As strongly emphasized by Charalam- 
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bous et al. , the obvious disadvantage of this technique is that ion bubbles "can- 
not provide the kind of detailed information about turbulent flow patterns that is 
available in the case of classical turbulence". We fully endorse this statement. Ac- 
tual visualization of quantized vortex filaments is also hardly possible (but see 
some new developments discussed below). On the contrary, solid, micron-size 
particles being used in the PIV and particle tracking techniques, albeit disturbing 
the flow, can visualize the normal flow and can even be used to "paint" quantized 
vortices. However, as was demonstrated in Sees. I2.2TI and 1275.31 in the low tem- 
perature regime, as T — > solid particles do not follow the flow, and, as was shown 
in Sec. 13.31 most likely cannot be trapped on quantized vortices. 

In the remainder of this chapter we will briefly discuss two new promising 
techniques that potentially can be used for direct flow visualization. 

7.2 Cavitating electron bubbles and metastable He2 molecules 

The first of these methods is being developed at Brown University by Maris and 
his colleagues 97 . The idea of the method comes from the fact that an application 
of a negative pressure to 4 He leads to an increase of the radius of the ion bub- 
ble which explodes, and its size grows substantially when the negative pressure 
has reached the critical value, P c . Furthermore, Maris et al. argue that since the 
pressure around the quantized vortex is reduced due to the Bernoulli effect, the 
ion bubbles trapped on the vortex have the size larger than those in the bulk of 
4 He and explode at the critical pressure of smaller magnitude (hence the possi- 
bility of direct distinguishing the bubbles trapped on vortices from those in the 
bulk of helium). In these experiments the pressure variations were produced by 
an acoustic, ultrasonic transducer generating either focused or planar (as in the 
last of cited works) sound wave. By choosing a suitable frequency and amplitude 
of the emitted sound, bubbles can be make so large that they can be visualized 
by the conventional optical methods (e.g. by a laser beam and a photomultiplier). 
The other advantage of this method is that it uses neither emission of the electron 
beam nor the electric field to accelerate the ions. The charged particles are those 
that already present in the bulk of 4 He, either as a result of direct ionization by 
e.g. cosmic rays, or through the more complicated process involving ionization, 
production of UV photons, and ejection of electrons from the walls into helium 
by the photoelectric effect. 

The absence of the applied electric field means that electron bubbles, unless 
they are trapped on quantized vortices, will faithfully follow the normal flow. In- 
deed, the results of visualization showed that in the weak, laminar counterflow 
most of the bubbles moved along the streamlines of the normal fluid, but along- 
side those about 10% of the observed electron bubbles followed snakelike paths; 
these could be assumed the electron bubbles trapped by vortices and sliding along 
vortex filaments. 

So far no new information on quantum turbulence has been obtained using this 
technique, but its relative simplicity has a potential to make it a tool competing 
with the PIV and particle tracking techniques. 

Another, completely new visualization technique is currently being developed 
by McKinsey and his co-workers 98 at Yale University. This technique employs 
metastable He2 triplet molecules, produced e.g. by the radioactive source, which 
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can be excited by two infrared photons from the ground a 3 E„ state to the (Pl,^ . 
About 90% of the excited molecules will decay into b^TIg state emitting photons 
at 640 nm well separated in the wavelength from the excitation photons and hence 
can be detected by standard techniques. All excited molecules decay back to the 
ground state within about 50 ns, so that the process can be repeated many times to 
make possible a detection of a single molecule. 

The authors claim that He2 molecules should be unaffected by quantized vor- 
tices at temperatures above 1 K and so allow the resolution of the normal flow 
even at the Kolmogorov scale. However, Vinen commented— that at sufficiently 
low temperatures triplet molecules may as well become trapped on quantized vor- 
tex core. In this case three-dimensional images of vortex lines can be produced by 
e.g. stereoscopic imaging. 

If successful, the further development of this method may be useful for recov- 
ering the velocity field and probability distribution function of the normal velocity 
fluctuations, and can also provide a new, nonintrusive tool for visualization of 
quantized vortices. 

8 Conclusions 

This article reviewed the recent progress in understanding the motion of solid par- 
ticles and mechanisms of their interactions with both the normal fluid and quan- 
tized vortices in turbulent 4 He. This problem is addressed in the context of the 
PIV and the particle tracking techniques recently implemented in 4 He, although 
we were certainly biased (perhaps, not surprisingly, considered our own research 
expertise) towards theoretical and numerical works and interpretation of experi- 
mental results. 

The problem of the interaction of particle tracers with the normal fluid and with 
quantized vortices is difficult, at least in principle: it is three-dimensional, time 
dependent and strongly nonlinear. This is why the best model available (the two- 
way coupling model described in Section 3) is computationally very expensive. In 
this review we have shown that, fortunately, some aspects of this particle-vortex 
interaction can be understood using simpler models, either numerical or analytical, 
particularly given the insight gained from the two-way coupling model. 

Although implemented recently (the first publications date 2002), the PIV and 
the particle tracking techniques applied to the problem of quantum turbulence in 
4 He had already produced some new results. Among them: the discovery of large- 
scale eddies, which do not have an analogue in the classical fluid dynamics, in the 
thermal counterflow past a cylinder; the wider and flatter, compared to classical, 
velocity profile in the turbulent boundary layer; scaling for the time evolution of 
the parameter characterizing reconnections of quantized vortices; non-Gaussian 
probability distribution function, with the tail scaled as v -3 , of the velocity i de- 
caying quantum turbulence. These results seem to suggest that both the PIV and 
the particle tracking techniques have great potential for studying quantum turbu- 
lence. Besides, the motion of solid particles in turbulent 4 He seems to be itself an 
interesting and non-trivial phenomenon worth the detailed theoretical and experi- 
mental study. 

Although the review of this rapidly development area of research cannot be 
comprehensive, we attempted to resolve, where possible, contradictions in theo- 



45 



retical interpretations of recent experimental observations and to address yet un- 
resolved issues. 

We also reviewed the related experimental methods, first of all the ion trap- 
ping technique which is perhaps the main tool for studying quantized vortices in 
the low temperature limit. We argue that these two experimental methods, i.e. the 
PIV/particle tracking and the ion trapping techniques do not compete but com- 
plement each other. Indeed, each of these techniques is most efficient in its own 
temperature and flow regime. 
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